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ABSTRACT 



Context. PIONIER stands for Precision Integrated-Optics Near-infrared Imaging ExpeRiment. It combines four 1.8 m Auxilliary 
Telescopes or four 8 m Unit Telescopes of the Very Large Telescope Interferometer (ESO, Chile) using an integrated optics combiner. 
The instrument has been integrated at IPAG starting in December 2009 and commissioned at the Paranal Observatory in October 
2010. It provides scientific observations since November 2010. 

Aims. In this paper, we detail the instrumental concept, we describe the standard operational modes and the data reduction strategy. 
We present the typical performance and discuss how to improve them. 

Methods. This paper is based on laboratory data obtained during the integrations at IPAG, as well as on- sky data gathered during the 
commissioning at VLTI. We illustrate the imaging capability of PIONIER on the binaries 6 Sco and HIP1 1231. 
Results. PIONIER provides 6 visibilities and 3 independent closure phases in the H band, either in a broadband mode or with a 
low spectral dispersion (R=40), using natural light (i.e. unpolarized). The limiting magnitude is //mag = 7 in dispersed mode under 
median atmospheric conditions (seeing < 1", r > 3 ms) with the 1.8 m Auxiliary Telescopes. We demonstrate a precision of 0.5 deg 
on the closure phases. The precision on the calibrated visibilities ranges from 3 to 15% depending on the atmospheric conditions. 
Conclusions. PIONIER has been installed and successfully tested as a visitor instrument for the VLTI. It permits high angular 
resolution imaging studies at an unprecedented level of sensitivity. The successful combination of the four 8 m Unit Telescopes in 
March 201 1 demonstrates that VLTI is ready for 4-telescope operation. 

Key words. Techniques: high angular resolution - Techniques: interferometric - Instrumentation: high angular resolution - 
Instrumentation: interferometers 



1. Introduction 



The Very Large Telescope Interferometer (VLTI, Haguenauer 
|et al.|2010| ) is equipped with four Unit Telescopes of 8 m (UTs) 
and four relocatable Auxiliary Telescopes of 1.8 m (ATs). It of- 
fers a unique combination of interferometric imaging capability 
and high sensitivity. A nice review of the scientific opportunities 
is being published by Berger et al. (2011, sub. to A&A Rev). It 
emphasizes the need to recombine a large number of telescopes 
to provide reliable snapshot imaging capabilities. However, the 
current instrumentation suite can only handle two or three tele- 
scopes simultaneously. The next generation of facility instru- 
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ments that will combine four telescopes is not expected to be 
operational prior to 2014. 

Therefore, IPAGQ and its partners have proposed to use a 
visitor instrumen^hdii provides VLTI, since 2010, with a new 
observational capability combining imaging, sensitivity and pre- 
cision. The principle of this Precision Integrated-Optics Near- 
infrared Imaging ExpeRiment (PIONIER) was approved by the 
ESO Science and Technical Commitee in spring 2009. The in- 
strument has been integrated at IPAG starting in December 2009 
and commissioned at the Paranal Observatory in October 2010. 
Only a few nights of commissioning were necessary before the 
instrument started routinely delivering scientific data with the 
ATs. On 17 March 2011 light collected by all four of the 8 m 
UTs was successfully combined for the first time. A picture of 
PIONIER installed at the VLTI laboratory is shown in Fig. [I] 
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Fig. 1. Picture of the PIONIER instrument installed in the VLTI 
laboratory in October 2010. The different modules refer to the 
conceptual description of Section [3] and Fig. [2] The size of the 
optical table is approximately 1.5 x 3 m. 



PIONIER relies on the integrated optics (10) technology to 
combine four beam s ( [Maibet et aL]|1999[ |Berger et aL]|200H 
Beni sty et al.|2009| ). In one observation, it provides the measure- 
ments of six visibilities and three independent closure phases 
with low spectral resolution (R « 40) in the H band. In many 
respect, it is complementary to the current AMBER instrument 
that provides three visibilities and one closure phase with vari- 
ous sp ectral resolutions ( R « 30, 1 500, 12 000) in the H and K 
bands ( [Petrov et al.|2007j ). PIONIER paves the road for the forth- 
coming instruments GRA VITY and MATISSE ( [Gillessen etaL] 
2010| |Lopez et al.|2008| ), by enabling the 4-telescope operation 
of VLTI. 

The paper is organized as follow: section [2] presents the 
main astrophysical programs of PIONIER and the correspond- 
ing high-level specifications. Section [3] describes the instrument 
and its different subsystems. Sectionj4] presents the PIONIER 
operations, including regular day-time maintenance, flux injec- 
tion, scientific fringe tracking and regular calibrations. Section|5] 
presents the data reduction strategy, focusing on the aspects that 
are specific to the instrument, or that have not been published 
elsewhere, and the current performance. Finally, Section [^illus- 
trates the imaging capability of PIONIER on the Be star S Sco 
and the spectroscopic binary HIP 11 231. The paper ends with 
brief conclusions and remarks about possible evolutions of the 
instrument. 



2. Science drivers and specifications 

PIONIER is designed for aperture synthesis imaging with a spe- 
cific emphasis on fast fringe recording to allow precision closure 
phases and visibilities to be measured. Among all the possible 
science cases, we defined three key astrophysical programs. 



2.1. Young Stellar Objects 

Young stars are surrounded by an accretion disk made of gas and 
dust, thought to be the birthplace of planets. Studying the inner- 
most astronomical unit is of major importance to understand the 
physical conditions for star and planet formation. 

Herbig Ae/Be stars are young stars of intermediate-mass, 
whose disks have been intensively studied with near-infrared in- 
terferometers in the last 13 years. Measurements enabled the de- 
termination of near-infrared sizes, found not to agree with clas- 



sical a ccretion disk model predictions (see Millan-Gabet et al. 
(2007 ) for a review). Consequently the view on the inner struc- 
ture of accretion disks has been substantially revise d. The con- 
cept of "puffed-up" inner r im has been introduced (Dullemond 



[etaL][20CTT| jlsella & N atta 2005). It marks the boundary in the 



disk where dust can no longer exist due to sublimation and 
forms a barrier to the photons potentially causing an alteration 
of the vertical structure of the disk. Recently, strong evidence 
was found that the matter flowing through this limit to the star 



is a major contributor to t he excess emission (Tan nirkulam et al. 
2008 ; Benist y et al.|2010 ). Revealing the nature of the matter in- 



side the dust sublimation limit and the exact structure of the in- 
ner disk at this transition is now investigated with care. We have 
identified a small sample of H Ae/Be stars whose environments 
have been clearly resolved by the VLTI and should therefore 
be ideal candidates for imaging with PIONIER in order to re- 
veal the morphology of the inner astronomical unit near-infrared 
emission. The simultaneous combination of four telescopes is re- 
quired to obtain reliable snapshot i mage reconstruction s of these 
complex and time-variable objects {Benisty et al. 2011). 



In addition, PIONIER will probe the structure of disks 
around the T Tauri stars, which are the low-mass analogs of 
Herbig stars. In these fainter objects, the inner edge of the dusty 
disk is closer to the star, and is not expected to be fully resolved, 
preventing us from direct imaging. Howe ver, the visibilities will 
give access to the inner disk location (Ake son et al.|2 005 



et al.||2009] ) and the closure phases will probe the morphology 



Eisner 



and surface brightness profile. In the best cases, it will be pos- 
sible to constrain, through proper radiative transfer modeling, 
the inner rim position, ellipticity and thickness. All of these are 
of major importance to understand the conditions for star and 
planet formation. As T Tauri stars are the faintest objects of 
our key programs, they set the need for limiting magnitude of 
//mag > 7. 



2.2. Debris disks around main sequence stars 



The inner solar system contains a cloud of small dust grains cre- 
ated when small bodies -asteroids and comets- collide and out- 
gas. This dust cloud, called the zodiacal disk, has long been sus- 
pected to have extrasolar analogs around main sequence stars. 
The very inner part of these exozodiacal disks have remained 
elusive until recently, when near-infrared interferometric obser- 
vations revealed a resolved excess flux of about 1 % around sev- 
eral nearby stars ([Absil et al.||2006| [2008] [20091 [oTFolco et al.| 
2007 j | Akeson et al .|2009| ). These excesses have been interpreted 



as being due to hot, possibly transient exozodiacal dust within 
1 AU from these stars. 

Our main goal is to survey main sequence stars to study the 
hot dust content during the late phases of terrestrial planet forma- 
tion, and its connection with outer resevoirs of cold dust. We ex- 
pect that about 100 main sequence stars of various ages between 
lOMyr and a few Gyr will be observed during the PIONIER 
lifetime, providing us with a solid statistical sample to derive the 
time dependence of the bright exozodi phenomenon. In addition, 
for the brightest cases, PIONIER could be used to characterize 
the grain properties thanks to multi-colour information, and to 
study the inner disk morphology. The challenge of this program 
is to obtain an accuracy of 1% or less on the calibrated visibility. 
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Table 1. Summary of the scientific requirements compared to 
the already demonstrated performance. Error is defined as the 
accuracy for one calibrated data point. 



Topic 


Sp. Band 


Mag. 


V 2 error 


CP error 


Herbig AeBe disks 


H,K 


> 5 


5% 


5deg 


T Tauri disks 


H,K 


> 7 


5% 


2deg 


Debris disks 


H,K 




1% 


ldeg 


Faint companions 


H,K 


> 6 




0.5deg 


Hot Jupiters 


H,K 






0.05deg 


Demonstrated 


H 


7.5 (AT) 


15-3% 


0.5deg 



2.3. Binaries, faint companions and Hot Jupiters 

Long baseline interferometric observations naturally comple- 
ment the radial velocity and adaptive optics surveys to discover 
and characterize binary systems. For instance, radial velocity 
measurements on massive stars are challenged by the lack of 
spectral lines and their intrinsic broadening. On active stars, the 
intrinsic radial velocity jitter may easily hide the signal of a 
faint companion. Moreover, compared with adaptive optics, the 
access to smaller separations and therefore shorter periodes in- 
creases the possibility to determine the dynamical masses of the 
components. 

Our first attempts with PIONIER (Absil et al. 2011, sub- 
mitted) demonstrate that the angular separation range from 5 to 
100 mas can be covered within one hour of observation, thanks 
to the simultaneous measurement of 4 closure phases at a low 
spectral resolution. Accordingly, we have initiated several sur- 
veys from massive to low mass stars, including young stars and 
stars in nearby moving groups. In order to reach a sufficient 
number of targets, these programs typically require a limiting 
magnitude of Hmag = 6 and a precision 6>/0.5deg on the clo- 
sure phases (1:200 dynamic range). 

Deep integration on few selected targets will permit to as- 
sess what is the best possible dynamic range achievable, until 
possibly reaching the planetary regime (1:2000). If successful 
PIONIER will have the means to separate the flux contribution 
of the planet from that of its parent star and potentially obtain 
low spectral resolution information. In order to achieve such a 
result, systematic biases on the closure phase estimate have to 
be hunted down and understood at the unprecedented level of 
0.05 deg ftZhao et al.\20W]\20Tl) . 

2.4. High level specifications 

The specifications of PIONIER have been defined as a trade- 
off between the requirements for the key astrophy sical programs 
and the technical contraints due to the project small size and 
fast timeline. A summary is presented in Table [T| We used the 
experience gained with ION IC-3 ( [Berger et al.]2003] ), VINCI ([Le] 
|Bouquin et al.|2004| 12006] ), and AMBER to make the following 
strategic choices: 

- the combination of 4 telescopes (6 baselines) simultaneously, 

- the use of an integrated optics beam combiner, 

- operation in the H- and K-band (not simultaneously), 

- a small spectral dispersion across few channels, 

- a fast readout of the camera. 

As built, PIONIER matches all these requirements except for 
the K band extension which is contemplated in the future. The 
instrument has been integrated and commissioned with the H- 
band integrated b eam combiner that w as already available at the 
IPAG laboratory ( [Benisty et al.|2009| ). 



3. Instrument description 

Figure [2] summarizes the key elements of the PIONIER instru- 
ment. Functionally, it consists of the following subsystems. 



3.1. Injection and Optical Delay Unit (IOPDU) 

The Injection and Optical Delay Unit (IOPDU) injects the free- 
space beams from the VLTI into optical fibers. It includes a tip- 
tilt correction, an OPD modulation and a polarization control. 
The IOPDU is made of four strictly identical arms, one per VLTI 
beam. 

1 . A dichroic mirror extracts the required band from the VLTI 
beams, while other wavelengths are used to feed the infrared 
VLTI guiding camera IRIS. The translation stage that sup- 
ports these optics has two observing positions (the H-band 
and K-band dichroic) and a free position to let the beam pass 
unaffected when PIONIER is not in use. 

2. A modulation of the optical path length is introduced by a 
mirror mount in a Physik Instrument (PI) piezo translation 
stage. It provides an OPD range of ±400 jim. 

3. The tip-tilt mirror, mounted on PI piezo devices, allows the 
beam angle to be corrected up to a frequency of 100 Hz with 
a range of ±200" (laboratory angle, corresponding to ±13 
times the size of the magnified PSF). 

4. A Lithium-Niobate plate of 2 mm thickness is used to com- 
pensate for the polarization phase- shift arising in the single- 
mode fiber of the IOBC (see Sec. |3.2| ). The amount of phase- 
shift between the vertical and horizontal axes is accurately 
adjusted by tilting differentially the plates in the individual 
beams. 

5. The off-axis parabola focuses the light into the single mode, 
polarization maintaining fiber. 

PIONIER' s fast tip-tilt correction comes in addition to the 
one provided by the telescopes. It mainly compensates for the 
additional contributions coming from the tunnel turbulence. The 
IRIS guiding camera provides the beam angle of arrival mea- 
surements at a rate from 100 H z to 1 Hz depending on the target 
brightness ( [Gitton et al.|2004| ). 

3.2. The Integrated Optics Beam Combiner (IOBC) 

The Integrated Optics Beam Combiner (IOBC) takes as input the 
signals from the IOPDU, and delivers 24 interferometric outputs. 

The IOBC is fed by polarization maintaining single-mode 
H-band fibers, whose lengths have been equalized with an ac- 
curacy of 20yum. Consequently, these fibers introduce a negli- 
gible differential chromatic dispersion. The polarization phase- 
shift between the vertical and horizontal axes (i.e. the neutral 
axes of the IO chip, aligned with those of the fibers) is of the 
order of 1 fringe. This polarization phase- shift can be compen- 
sated by adjusting the inclination of the Lithium-Niobate plate 
in the IOPDU (see Sec.[3J}. 



The design of the IOBC is described in detail in Benisty 
et al.| ( 2009| ). The four incoming beams are split in three and 
distributed in the circuit in a pairwise combination. A so-called 
"static-ABCD" combining cell is implemented for each base- 
line. It generates simultaneously four phase states (almost in 
quadrature). Consequently, 24 outputs have to be read. The low- 
chromaticity phase shift is obtained through the use of specific 
waveguides with carefully controlled refraction index. This com- 
biner can be used both in fringe- scanning mode (VINCI-like) or 
AB CD-like mode. There is currently one combiner available in 
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Fig. 2. PIONIER conceptual scheme. The Injection and Optical Delay Unit (IOPDU) is responsible for the injection of the four 
beams (either from CAU or from VLTI) in the single mode fibers. It also implements the OPD modulation, the tip-tilt correction and 
the polarization control. The Integrated Optics Beam Combiner (IOBC) recombines the four beams into 24 interferometric outputs 
in a pairwise scheme. These outputs are then imaged, and eventually spectrally dispersed, by the Imaging Optics and Dispersion 
Unit (IMODU) on the detector (DET). The aim of the Calibration and Alignment Unit (CAU) is to provide four interferometric 
beams that match the size and angle of the VLTI beams. 



the H ban d and one that is being developed for the K band ( Jocou 
|et al.|2010) . 

3.3. Imaging Optics and Dispersion Unit (IMODU) 

The Imaging Optics and Dispersion Unit (IMODU) images the 
24 outputs of the IOBC onto the camera's focal plane, with or 
without spectral dispersion in the perpendicular direction (1, 3 or 
7 spectral channels across the H band). This unit is the one that 
was used with the IONIC- 3 instrument at IOTA (Ber ger et al. 
2003). Its image quality permits to focus about 80% of the flux 
into a single pixel of the PICNIC camera (pixel size 40yum). A 
Wollaston prism can be inserted to acquire separately the fringes 
for the two linear polarizations (vertical and horizontal axes). 
This device is only used to internally adjust the compensation 
of the polarization phase- shift (see Sec. |3.2| ). Once this adjust- 
ment is done, the Wollaston prism is removed and the fringes are 
acquired in natural (unpolarized) light without loss of contrast. 
Consequently the instrument provides no polarization informa- 
tion in targets. 



In practice, these prisms are inserted and removed manually 
from the optical beam. The positioning of the Wollaston prism is 
repeatable at ±0.1 pixel. The positioning of the dispersive prism 
is repeatable at ±0.1 pixel. Altogether, these measurements are 
compatible with a stability of ±1% for the effective wavelength 
and for the flux splitting ratio between the outputs. 



3.4. The detector (DET/DFE) 

The detector (DET) comprises the focal plane array, the cryostat, 
the filter wheels, and the internal electronics. More electronic 
functions, including the digitizing of the video signal, are pro- 
vided by the Detector Frontend Electronics (DFE). The d etector 
is the PICNIC camera previously used at PTI and IOTA (Traub 
et al. ||2004} |Pedretti et al.||2004| ). Its electronics has been par- 
tially upgraded to improve the frame rate. The array is clocked at 
5 MHz and performs analog-to-digital conversions every 0.4 jus. 
The main limitation for speed is the time for the video signal to 
stabilize once a pixel has been addressed 3 //s). The resulting 
time for the elementary operations reset and read are presented 
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Table 2. Typical time to execute elementary operations with our 
fastest possible clocking of the PICNIC camera. These numbers 
are for a single analog-to-digital conversion per pixel. 



Detector elementary 


Required time 


operation 


Qis) 


Reset full quadrant 


205 


Read 24 pixels x 1 line 


190 


Read 24 pixels x 7 lines 


1150 



in Table [2] Consequently, the detector achieves a frame-rate of 
0.85 kHz when reading 168 pixels (24 outputs times 7 spectral 
channels) in a non-destructive way. 

We characterized t he readout nois e follo wing the same pro- 
cedure as described by |Pedretti et al.|p004| ). We found a noise 
of « I3e~ per read while Pedretti et al. reported « 9e~ per read 
on this camera. This noise is reduced to « 8e~ by averaging 8 
successive analog-to-digital conversions of the same pixel. It is 
still possible to reduce the noise by averaging several frames in 
a non-destructive way, although it has a strong impact on the 
minimum integration time possible. 

3.5. The Calibration and Alignment Unit (CAU) 

The Calibration and Alignment Unit (CAU) serves dual func- 
tions: (1) to inject, as a substitute for the astronomical beams, 
mutually coherent beams derived from a common Tungsten 
source, for laboratory verification and health check purposes; (2) 
to reflect, via a corner cube mirror, a reverse-propagating beam 
from the output of the IOBC towards the IRIS guiding camera, 
for the purpose of verifying the alignment of the PIONIER in- 
strument with the IRIS reference positions. 

3.6. Control System and electronics (CS) 

Control System (CS) includes hardware control of the instru- 
ment units, detector readout, quicklook and interaction with 
VLTI. The hardware has been designed to follow the standard 
VLTI architecture including an instrument workstation (WS), an 
Instrument Control System (ICS), a Local Control Unit (LCU) 
and a Detector Control System (DCS) running on the detector 
WS. The instrument WS is a standard Linux PC platform. The 
detector WS is an industrial Linux PC platform equipped with a 
PCI-7300A board used to generate clock signals and to acquire 
detector data. Both stations are running under the VLT Software. 

Instead of the traditional LCUs, the instrument is controlled 
by an Embedded PC and several EtherCAT modules like remote 
10 and stepper motion controllers from the company Beckhoff. 
Through a fruitful collaboration with ESO, PIONIER has been 
used as a pilot project for testing these commercial components 
and developing the associated VLT Software extensions. This 
success not only tries to simplify the coding of ICS but it also 
attemp ts to introduce tec hnologies envisaged for E-ELT instru- 
ments ( Kiekebusch et al. 2010 ). 

PIONIER electronics is housed in a cooled cabinet close to 
the instrument optical table. It includes the detector front-end 
electronics, the DC detector power supplies together with the 
detector WS, tip-tilt and scanning piezo controllers, EtherCAT, 
calibration lamps and control electronics for electro-mechanical 
devices. 



3.7. Interfaces 

Dialog between PIONIER and VLTI (e.g. sending new coordi- 
nates) is achieved through the VLTI Interferometer Supervisor 
Software. The detector WS is also listening to the VLTI 
Reflective Memory Network (RMN) to get real-time centroid 
information from the IRIS guiding camera. 

The raw output data follow the standard defined by ESO for 
its interferometric instruments (Ballester & Sabet 2002), includ- 
ing the general FITS header. The only difference is the addition 
of a new dimension (depth) in the column DATAj of the binary 
table IMAGING _DATA. The third dimension of DATAj repre- 
sents the consecutive analog-to-digital conversions of each pixel. 
Consequently, a typical night provides several gigabytes of data. 

The Observation Software running on the instrument work- 
station coordinates the execution of an exposure for a given 
observing mode. From the high level point of view, PIONIER 
is operated via the Broker of Observing Blocks that exe- 
cutes Observing Blocks (OBs) fetched from the standard p2pp 
ESO software. OBs can be conveniently generated by the 
aspro2^] preparation software from the Jean-Marie Mariotti 
Center (JMMC). 



4. Operating PIONIER 

This section describes the PIONIER operations, including reg- 
ular day-time optical alignment, flux injection procedure, scien- 
tific fringe recording and regular calibrations. 

4.1. Optical alignment 

As a visitor instrument, PIONIER is maintained and operated by 
the PIONIER team. It is only operated in visitor-mode and the 
observer is responsible for checking the health of the instrument 
at the beginning and during the observing run. The day-to-day 
stability of the beam angle inside PIONIER is « 4", that is a 
quarter of the magnified PSF size. This is well within the range 
of the tip-tilt actuators (±200"). The day-to-day position stabil- 
ity of the fringes acquired with the internal calibration is better 
than 15//m in OPD. The same typical variations are observed 
with the VLTI calibration source MARCEL. Again, this is well 
within the range of the modulation mirror (±500 jim) and does 
not impact the operations. 

However, the two following adjustments should be carefully 
checked at the beginning of each observing run: 

- the vertical/horizontal positioning of the IOBC, in order to 
optimize the coupling of the output spots into the detector 
pixels ; 

- the tilt angle of the Lithium-Niobate plates, as the polarisa- 
tion phase- shifts of the fibers drift by typically Xj 8 over a 
month timescale. 

These adjustments remain valid for the few following days 
of a typical run. 

4.2. Flux injection 

The stellar images are stabilized in tip-t ilt by the IRIS gu iding 
camera located in the VLTI laboratory ( [Gitton et al.|2004| ), only 
a few meters away from PIONIER. The image shift caused by 
the difference of atmospheric refraction index between the H- 
band (PIONIER fibers) and the K-band (IRIS guiding camera) 



3 http :// w w w. j mmc . f r/aspro _p age . htm 
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Fig. 3. Scanned interferograms obtained with PIONIER on the high SNR regime (left) and the low SNR regime (right) and using 
the spectral dispersion of 7 channels across the H-band. An example with the baseline with the smallest frequency is represented 
here. The bottom, black signal represents one output of the combiner, for the central spectral channel (i.e. one pixel of the detector), 
versus the optical path difference in the scan. The large temporal coherence and the photometric fluctuations are clearly visible. The 
top, red signal represents the subtraction of two conjugate outputs, after averaging the 7 spectral channels. The SNR is enhanced 
and the photometric fluctuations are almost perfectly removed. 



is handled by the Interferometer Supervisor Software. The soft- 
ware automatically offsets the IRIS guiding point so that the stel- 
lar image in the H-band is kept fixed whatever the airmass and 
pointing direction. 

The amount of flux injected in the fibers is optimized by a 
grid search with the internal tip-tilt of PIONIER. We use a square 
grid of 10x10 steps of 2.5" each. The resulting image is fitted by 
a 2D gaussian. The procedure is applied sequentially on the four 
beams, during the evening twilight, using the beacon sources lo- 
cated on the VLTI telescopes. The resulting accuracy is better 
than l/50 th of the PSF size. The precision is dramatically re- 
duced by speckle noise when the procedure has to be performed 
on stellar light. 



Table 3. Fringe frequencies versus the number of spectral chan- 
nels to be read and the number of steps in scan. The scan length 
corresponding to f m i n is set to 80 /mi. The modes are sorted from 
less sensitive (top) to most sensitive (bottom). 



Detector 


# of spectral 


steps in 


fmin 


steps per fringe 


mode 


channels 


scan 


(Hz) 


fmax 


Double 


1 


512 


160 


3.6 


Fowler 


7 


512 


80 


3.6 


Fowler 


7 


1024 


40 


7.2 


Fowler 


1 


1024 


240 


7.2 


Fowler 


1 


2048 


120 


14.4 



4.3. The interferometric signal 

The fringes are temporally scanned across the coherence length 
by means of the long-range OPD modulation devices. Figure [3] 
shows examples of raw interferograms obtained in the high SNR 
and low SNR regimes. The four beams are modulated at non 
redundant velocities (-3v, -v, +v, +3v) in order to generate the 
six fringe signals. The speed v is defined by the ratio between the 
length of the smallest scan (set to 80 jim) and the duration of a 
scan. For faint stars, its optimal value is a compromise between 
sensitivity and the necessity to freeze the atmospheric effects for 
the fringes at the lowest frequency f m { n - 2v/A. For bright stars, 
v is imposed by the detector frame rate and the need to correctly 
sample the fringe at the highest frequency f max = 6v/A. 

The most used scan setups are summarized in Table [3] In 
mode "fowler", the detector is reset once and is left integrating 
non destructively during the entire scan. This mode achieves the 
best frame rate and the lowest readout noise (when oversampling 
the fringes). In mode "double", the detector performs a reset- 
read-read sequence for each step in scan. This mode is used on 
bright stars that would saturate the detector in "fowler" mode. 
Its main drawback is to slow down dramatically the frame rate. 
Consequently, it is not possible to use the mode "double" with 
the spectral dispersion over 7 channels as the resulting frame rate 
would be too slow compared with the atmospheric coherence 
time. 



4.4. Quick-look and fringe tracking 

During the observation, each scan is processed by a quick-look 
algorithm that also implements a slow group-tracking technique: 

1 . The spectral channels are added together in a pseudo broad- 
band signal. The four outputs of each baseline are summed 
together taking into account their relative phase shift. The re- 
sult is a single interferometric signal per baseline represented 
in red in Fig. [3] 

2. At the end of the scan, the signal is processed by the algo- 



2005). It provides a SNR 



rithm presented by Ped retti et al. 
and a group-delay position for each baseline. 

3. Then these six group-delay measurements are converted into 
offsets for the four input beams. The inversion algorithm 
only takes into account the baselines with an SNR larger than 
3. The redundancy is used to recover the offsets of a maxi- 
mal number of beams. The results are an optimal estimate of 
the offsets to be applied to the beams and a detection flag per 
beam (0/1). 

4. For all beams with a positive detection flag, the zero posi- 
tions of the scanning piezo-electric devices is updated with 
the new offset value, before the next scan starts. 

The typical repetition rate of the loop (1 Hz) is enough to 
keep the group-delay residuals within ±20 yum under median at- 
mospheric conditions. The redundancy of the 6 baselines to re- 
cover 3 optical path differences makes the tracking robust against 
temporary flux losses or strongly resolved baselines. 



6 



J.-B. Le Bouquin et al.: PIONIER: a 4-telescope visitor instrument at VLTI 



4.5. Observations and calibration 



5.1. Cosmetic 



Each pointing typically produces 5 files of 100 fringe scans. 
Then a dark exposure is recorded in a file of 100 scans with the 
same detector setup but with all shutters closed. The flux split- 
ting ratios of each input are finally calibrated by files of 50 scans 
with one beam illuminated at a time. 

The sequence takes between 5 min and 10 min depending on 
the modulation speed. An additional 3 min are needed to pre- 
set the VLTI telescopes and delay lines and acquire the star. A 
single calibrated point (science + calibrator) takes between 15 
to 30 min. It is possible to skip the optimisation of the telescope 
guiding system if the V magnitudes of the science and calibra- 
tion stars are within 0.5 mag. It significantly speeds up the target 
re-acquisition (1 min instead of 3 min). 

The spectral calibration is performed at the end of the night. 
It consists of long scans of 200 yum with high SNR recorded 
in the internal CAU. The effective wavelength is estimated by 
Fourier Transform Spectrometry. We checked the linearity of the 
OPD modulation by observing a 1 .55 jim infrared laser. The day- 
to-day repeatability of the effective wavelength is ±1%, which 
is consistent with the optical stability presented in the previous 
section. Taking into account the systematic uncertainties, the fi- 
nal accuracy on the effective wavelength is ±2%. It matches the 
requirement of the key programs presented in Sec. [2] In particu- 
lar, it meets the required 1% accuracy on the calibrated visibility 
associated to the most demanding program, for which the visi- 
bility mo dels are not very sensitive to the wavelength calibration 
(Sec |2.2| faint debris disks around marginally resolved photo- 
spheres). 



5. Data reduction and performances 

In this section, we focus on the steps of the data reduction that 
are specific to the PIONIER instrument, or that have not already 
been published. 

The data reduction software of PIONIER converts the raw 
FITS file produced by the instrument into calibrated visibilities 
and closure phase measurements written in the standard OIFITS 
format ( [Pauls et aL]|2005| ). These output files are science-ready 
and can be directly hand led by software such as LITpro (model 
fitting, Tallon-Bosc et al. 2008) or MIRA (image reconstruction, 
Thiebaut 2008J). The PIONIER data reduction software is written 
as a yori ck^jpackage called pndrs and publicly availabl^] 

Our strategy is directly inspired from the FLUOR, VINCI 
and IONIC- 3 experiments. In these instruments, the scanning 
method associated with an estimate of the visibility in the 
Fourier space provided sensitive and accurate measurements, 
without the need for tricky internal calibrations. Basically, the 
intensity measured at the output A (respectively B, C and D) of 
the baseline ij of the combiner is written as: 



i l l=K l lP l +K ] lP ] + 



ylkVqPPjfi V J COS(-- + tft + 0")(D 



where ij are the baselines (12, 13, 14, 23, 24, 34), F are the 
photometric flux from input beams, \i 3 A and (p A are the instru- 
mental visibilities and phases of the beam combiner, k are the 
internal flux splitting ratio, and S lJ is the OPD modulation across 
the scan plus the atmospheric piston. The complex visibility of 
the observed target is V lJ e l ® 13 . 



4 http://yorick.sourceforge.net 

5 http ://app s j mmc.fr/~swmgr/pndrs 



The basic operations that are related to the PICNIC detec- 
tor are performed first. The consecutive non-destructive reads 
are subtracted accordingly to the detector mode ("double" or 
"fowler") to provide the intensity. The dark level is removed. 
In the "fowler" mode of the detector, a dark level is estimated 
independently for each step in scan. 

The flux splitting ratios are estimated from the dedicated cal- 
ibration files. They are computed independently for each wave- 
length. We use these flux splitting ratios to compute what we call 
the flat-fielded kappa matrix: 



and the flat-fielded signal: 



7*7 



(2) 



(3) 



We use the instrumental visibilities \i 3 A coming from laboratory 
measurements of the beam combiner, actually close to one. 



5.2. Disentangling the photometric and interferometric 
signals 

The pairwise design allows the instantaneous photometric signal 
to be extracted by a global fit to the coherent data. The main 
idea is to add together the outputs A,B,C and D from a given 
baseline, after flat-fielding: 

c ij = ^ + ?j + ^ + ?i (4) 

The opposite phase relation between A and C (respectively B and 
D) due to energy conservation removes the interferometric part 
of the signal: 



(5) 



n. A , a j t , ,v c , k'^. Equation [5] defines a linear sys- 
tem between the six c lJ signals (one per baseline) and the four P l 
photometric signals (one per beam). This system can be inverted 
and is well constrained thanks to the pairwise combination. We 
obtain an estimate of the P l for each step in scan. This method to 
recover the photometric sign als was used in the P IONIER pre- 
cursor instrument IONIC-3 ( [Monnier et al.||2004| ). Although it 
is not widely known, this method is a fundamental a dvantage of 
multi-telescope pairwise combination (see $4.1.2 in |B find et al 
|20TTT ). 

We build a single interferometric signal per baseline com- 
bining the information of all outputs. We first inject the P l esti- 
mates into the flat-fielded kappa matrix in order to remove the 
continuum part of the interferometric signals. Then we co-add 
the A,B,C and D outputs taking into account their relative phase 
shift (estimated from laboratory measurements): 



Z # 



(6) 



o=ABCD 



5.3. Closure phase estimate 

The interferometric signal m lJ is filtered in Fourier space in or- 
der to keep only the frequencies near the modulation frequency 
of each baseline. Then the baselines are combined together into 
bispectra, that are integrated over the length of the scan: 



opd 



(m ik T 



(7) 
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Fig. 4. Examples of statistic of the bispectrum over 100 scans ob- 
tained in the high SNR (left) and low SNR (right) regimes. The 
horizontal and vertical axes are the real and imagery parts of the 
bispectrum respectively (arbitrary unit). The phase of the bis- 
pectrum is free from fringe modulation and atmospheric piston 
terms, and points toward the closure phase. The closure phase 
average (solid line) and standard deviation (dashed lines) are es- 
timated with the bootstrapping technique. Figure [3] contains typ- 
ical interferograms corresponding to these measurements. 



Four bispectra are formed Z? 123 , b 124 , b 234 and b 134 . Figure [4] 
shows examples of bispectrum statistics obtained for 100 scans 
in the high SNR and low SNR cases. The phase terms related 
to the fringe modulation and the atmospheric piston cancel out 
and all bispectrum estimates point toward the same phase: the 
closure phase ® u + O 71 - <£> lk . We take the argument of the bis- 
pectrum average to estimate the final closure phase. We use a 
bootstraping method to compute its standard deviation. 



5.4. Visibility estimate 

The critical step is the photometric normalization of the interfer- 
ometric signal m lj by the photometric signal VP Z P J . We imple- 
ment two versions. 

For bright stars, the photometric normalization is performed 
independently for each step in scan. This simultaneity is the key 
to obtain prec isely calibrated interfero grams on bright stars. As 
discusse d by[Coude Du Fo resto et al.| ( |1997) an d |Kervella et al.| 
( 2004b ), we use a Wiener filtering of VP'PA Still this is not 
enough to avoid numerical instabilities when l pj get close to 
zero. We filter out the scans with such a behavior. 

For fainter stars, the photometric signal is averaged over the 
length of a scan. When dominated by the detector noise, it dra- 
matically enhances the precision on the photometric signal. In a 
sense this is an extreme case of the previous method, when the 
Wiener filtering is so strong that it is equivalent to a simple time 
average. 



After the photometric normalization, we use the method de- 
tained by Coude Du Foresto et al. (1997), applied independently 
to each spectral channel and baseline, to obtain an estimate of 
the square visibility |V U | 2 . This method is based on the Power 
Spectral Density (PSD) of the signal to disentangle the fringe 
power (peak) from the detector and photon noise biases (con- 
tinuum). Typical PSD are presented in Fig. [5] The shape of 
the underlying bias power can be correctly fitted with the form 
a + b sm 2 (2ncr/cr max ) where cr is the wavenumber. The first term 
a corresponds to the photon noise and the second term b cor- 
responds to the detector noise. This underlying structure of the 
noise is due to the "fowler" readout scheme: the detector noise 



is high-pass filtered by the difference between two consecutive 
reads. 



5.5. Transfer function calibration 

The instrumental and atmospheric contributions to the visibili- 
ties and closure phases (the so-called transfer function) are mon- 
itored by interleaving the observation of science targets with cal- 
ibration stars. We implemented several strategies to interpolate 
the transfer function at the observing time of the science targets. 
We found that the most robust is to simply average all the trans- 
fer function estimates of the night, with a weight inversely re- 
lated to the time separation. The uncertainty on the transfer func- 
tion at the time of the scientific observations takes into account 
both the error bars and the dispersion of the transfer function es- 
timates: a sequence of precise (small error bars) but dispersed 
points leads to large uncertainties on the calibrated product. An 
example of a transfer function with its interpolation is displayed 
in Fig.[p 

Calibration stars fo r bright targets are g enerally extr acted 
from t he catalogues of Meran d et al. (2005) and Borde et al. 
(2002). Calibration stars for f ainter targets are ext racted from 
the JSDC catalog from JMMC ( [Lafrasse et al.|2010 ). 



5.6. Performances 

In this section, the term median atmospheric conditions refers to 
seeing « 1", To ~ 3 ms, and wind speed « 6 m/s. 

The limiting magnitude of PIONIER is defined by its capa- 
bility to detect the fringes during the observations (SNR^3 per 
scan). The typical measured sensitivity is //mag ^ 5 for a dis- 
persion over 7 spectral channels and using a scan of 512 steps 
with the detector in mode "fowler". When the atmospheric con- 
ditions are better than the median, it is possible to reach a lim- 
iting magnitude of //mag « 7 using a scan of 1024 steps. In 
broad band, we have been able to detect and track the fringes on 
several unresolved stars with //mag > 8 under good conditions. 
These numbers are for the Auxiliary Telescopes. 

The statistical uncertainties on the closure phases (error bars 
in Fig. [6j right) are compatible with the data dispersion. The 
final accuracy of the calibrated closure phases generally ranges 
from 5 deg to 0.5 deg. These uncertainties not only depend on the 
target brightness but also on the turbulence strength. Sequences 
with stable closure phases, down to 0.1 deg, have been recorded 
on bright stars under atmospheric conditions better than median. 
At such a level, new possible biases should be studied such as 
the dependence of the instrumental closure phase on the spectral 
type. This has not been done yet. 

The measured accuracy on calibrated visibilities ranges from 
15 to 3% depending on the atmospheric conditions. On bright 
targets, the dispersion of the uncalibrated data is sometimes 
larger than the statistical uncertainties computed by pndrs 
(Fig.[6j left). It means that we are sometimes facing some non- 
stationary biases. Our prime suspects are the mechanical OPD 
perturbations that could be related to the wind strength, perhaps 
by shaking the telescope structure. Indeed, when the wind speed 
is higher than 6 m/s, the PSDs recorded by PIONIER are clearly 
affected by vibrations, as seen in Fig.[5]left. These perturbations 
are typically in the range 10-100 Hz, with a typical life-time of a 
few seconds. This should be further investigated before drawing 
definite conclusions. 

The VINCI and IONIC- 3 instruments achieved a visibil- 
ity accuracy of 1%, using a similar instrumental concept to 



8 



J.-B. Le Bouquin et al.: PIONIER: a 4-telescope visitor instrument at VLTI 



o 

GO 



o 





Fig. 5. Examples of average Power Spectral Density of the interferometric signals (100 scans) obtained in the high SNR (left) and 
low SNR (right) regimes. Only a single spectral channel of one baseline is represented. The fringe signal corresponds to the peak at 
o~ = 0.6 /mi -1 (A = 1.66 /mi). The dashed red line is the fit of the underlying bias power by its theoretical form a + b sm 2 (2ncr/cr max ). 
In the low SNR regime, electronic artefacts can be seen at cr « 0.22 yu m" 1 , 0.8 yum -1 and 0.97 yum -1 . In the high SNR regime, the 
fringe peak is clearly aliased by mechanical vibrations in VLTI (see Sec. 5.6 ). Figure[3]contains typical interferogram corresponding 
to these measurements. 
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Fig. 6. Examples of transfer function calibration for the visibility (left) and the closure phase (right). Only a single spectral channel 
and a single baseline (triplet for th e closure phase ) is re presented. The black squares are the transfer function estimates made on 
calibration stars extracted from the Merand et al. ( |2005 ) catalogue. The solid and dashed lines represent the interpolation of the 
transfer function, with its uncertainty. The colored circles are the observations of scientific targets. 



PIONIER (see for instance |Kervella et^]|2004a| ). We plan to 
dedicate a few observing nights to define the best observing 
setup to reach the same level of accuracy. 

6. Image reconstruction of 8 Sco and HIP11231 

In this section we illustrate the imaging capability of PIONIER 
on the stars 6 Sco and HIP11231. S Sco is a well known bright 
Be star (//mag = 2.1) accompanied by a faint companion on a 
highly eccentric orbit. HIP11231 is a G1V spectroscopic binary 
with apparent magnitude //mag = 4.8. The PIONIER obser- 
vations are presented in Appendix [A] and [B] together with their 
(w, v) plane coverage. 6 Sco has been observed once with the 
intermediate configuration (D0-H0-G1-I1). HIP 11231 has been 
observed at 7 epochs, and with the three AT quadruplets (D0-H0- 
Gl-Il, E0-G0-H0-I1, A0-K0-G1-I1). Data have been processed 
with the pndrs package. 



6.1. Image reconstruction and model fitting of 5 Sco 



We perform model-in dependent image reconstruction with the 
MIRA software from Thiebaut (2008). The starting point is a 
delta function at (0, 0). The pixel scale is 0.2mas and the field-of- 
view is 200 x 200 pixels. In practice, the MIRA algorithm mini- 
mizes a joint criterion which is the sum of (1) a likelihood term 
which measures the compatibility with the data, and (2) a regu- 
larization term which imposes priors on the image. The relative 
weight between these two terms is controlled by a multiplicative 
factor called "hyperparameter". We use the "total variation" reg- 
ul arization associated w ith positivity constraint as recommended 
by |Renard et al. ( |2011| ). We set the hyperparameter with a small 
value of 100, so that the weight of the regularization term is kept 
small with respect to the fit to the data. It brings some superreso- 
lution, at the cost of an increased level for the noise in the image. 
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2011-06-04 



0.15 




— 0.05 



Fig. 7. Image reconstruction of 6 Sco. The reference frame is 
centered on the primary star. The red cross is the position of the 
com panion at the t ime o f the PIONIER observations predicted 
from Tycner et al. (201 1 ). Scale is in mas, North is up and East 
is to the left. The image maximum is normalized to unity and the 
color scale is linear from to 0.15. 

Table 4. Best fit disk+point model for 6 Sco: Uniform Disk (UD) 
diameter of the primary, flux ratio (p) and position (£ra and £dec) 
of the secondary. 



UD 


p 6m 


fidec xl 


[mas] 


[%] [mas] 


[mas] 


1.48 ±0.07 


6.3 + 0.5 12.72 + 0.1 


5.96 + 0.1 1.1 



We use the information of all the spectral channels to improve 
the (w, v) plane coverage ("gray image" hypothesis). 

The reconstructed image is presented in Fig. [7] The compan- 
ion is easily detected. The main component appears extended 
as the disk surrounding it is marginally resolved by the VLTI 
baselines. This result illustrates the snapshot imaging capability 
of PIONIER thanks to the simultaneous beam combination of 4 
telescopes. 

In parallel, we perform a fit with a binary model with one re- 
solved component. Results are presented in Table 4] The result- 
ing visibility and closure phase curves are overlaid to the data in 
Appendix [A] We checked that the position of the co mpanion is 
consistent with the prediction by Tycn er et"aL| ( |201l| ). This im- 
age reconstruction of S Sco demonstrates that the pndr s package 
computes the sign of the closure phases and of the (w, v) plane in 
a consistent manner. In other words, it means that the produced 
QIFITS files are c ompliant with the sky orientation defined by 
|Pauls etaL|j2005) . 

6.2. Image reconstruction, model fitting and orbital solution 
of HIP11231 

We use the same MIRA parameters as for 6 Sco except that the 
pixel scale is set to 0.5mas. The reconstruction is performed in- 
dependently for each epoch. We use the information of all the 
spectral channels to improve the (u, v) plane coverage ("gray im- 
age"). Reconstructed images are presented in Fig. [8] The spec- 



Table 5. Campbell orbital elements of HIP 11231 from our 
PIONIER astrometric measurements and from IJancart et al.l 
(20051. Both solutions adopt the spectroscopic period P = 



142.33d ( [Barker etal.|1967) . 



Parameters 


To 


a + 


e 


i 


a 


CO 




[d] 


[mas] 




[°] 


[°] 


[°] 


Jancart, 2005: 


37159.1 


8.09 


0.29 


60.6 


200.2 


188.2 


This work: 


37187.4 


21.3 


0.25 


57.5 


232.6 


194.8 



semi-major axis of secondary orbit in this work. 



troscopic binary is easily resolved by PIONIER+VLTI with the 
three AT configurations. The noise in the images is kept below 
2.5% of the image maximum. The flux ratio between the primary 
and the secondary, measured by aperture photometry in each re- 
constructed image, ranges from p = 0.18 top = 0.21. 

In parallel, we perform a global fit to the full data set with 
the following parameters: (1) the flux ratio, forced to be iden- 
tical for all epochs and (2) the relative positions, that we let 
free between each epoch. We combined together all the spec- 
tral channels to improve the (u, v) plane coverage. Such a global 
fit is conveniently performed with the LITpro software ( Tallon-| 
|Bosc et al.|[2008] ). The best-fit flux ratio is p = 0.206. At each 
epoch, the relative separation could be retrieved unambiguously. 
The fitted positions and the flux ratio are consistent with the re- 
sults from model-independent image reconstruction. The result- 
ing visibility and closure phase curves are overlaid to the data in 
Appendix |B 

shows the relative separations derived from the fit 
th the best orbital solution. We adopt the spectro- 
scopic period derived by |Barker et al.| ( |1967| ) from radial ve- 
locity, but we allowed all the other Campbell elements to vary, 
namely: the time passage through periastron Tq, the eccentric- 
ity e, the semi-major axis a, the position angle of the ascending 
node Q, the argument of periastron oj and the orbital inclination 
i. Parameters of the best orbital solution are presented in Table [5] 
togeth er with the solution determined from Hipparcos astrome- 
try by |Jancart et al.| ( [2005]). The di fference in semi-major axis a 
is due to the fact that ] Jancart et al.| (|2005) reconstructed the orbit 
of the barycenter of light from unresolved astrometry, while we 
reconstruct the orbit of the secondary from resolved astrometry. 
The difference in time passage through periastron To can be ex- 
plained by the uncertainties on the period and by the fact that our 
observations have been taken several orbital periods later than 
the Hipparcos ones. The disagreement in Q is more striking as 
we cannot find a simple explanation. 



Figure[9 
overlaid wit 



7. Conclusions and perspectives 

The 4-telescope PIONIER instrument has been successfully de- 
veloped and commissioned at the VLTI interferometer. The use 
of integrated optics technologies combined with the expertise 
gained from projects such as FLUOR, AMBER, VINCI and 
IONIC-3 made it possible to design, build and commission the 
instrument in about 18 months. PIONIER is routinely used for 
science operations in the H-band and permits high angular res- 
olution imaging studies at an unprecedented level of sensitivity 
and precision. The overall performances are in agreement with 
the requirements from the key programs, although some work 
is still needed to achieve the best possible accuracy on the cali- 
brated visibilities and closure phases. 
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Fig. 8. Model independent image reconstruction of HIP1 1231 with the MIRA software. Each night has been reconstructed indepen- 
dently. The reference frame is centered on the primary star, which saturates the color table. Scale is in mas, North is up and East is 
to the left. The image maximum is normalized to unity and the color scale is linear from to 0.15. 




da (mas) 

Fig. 9. Best fit of the apparent orbit of HIP11231, considering all the Campbell elements free except the period (142.33days, from 
Barke r et al.|1967j ). The black circles are the PIONIER astrometric measurements. The red dotted line is the best orbit with its line 
of nodes in dashed. The red crosses are the positions of the companion at the time of the PIONIER observations predicted by the 
best-fit orbit. The average residual distance is 60yuas. The reference frame is centered on the primary star. North is up and East is to 
the left. 



PIONIER is expected to stay a few years at VLTI as a visitor 
instrument, until the arrival of the second generation instruments 
MATISSE and GRAVITY. Within this lifetime, several improve- 
ments are contemplated: 

- The imaging and dispersion unit will be upgraded with 
a motorized translation stage for the dispersion prisms. It 
will allow to change the operational mode (broad band 
or dispersed) during the night. This should be installed in 
PIONIER around December 201 1. 

- Changing the operational wavelength from the H-band to the 
K-band would benefit all the key science programs. A com- 



ponent is already available at IPAG. The design of the as- 
sociated imaging optics is under study. Ideally, the changes 
will be installed in PIONIER early 2012. 
- The arrival of a new generation of infrared detectors should 
provide detector noise as low as 2.5 e~ at a frame-rate larger 
than 1 kHz in destructive mode ( |Gach et aL 2009 Finger 
|et al.||2010| ). Installed in PIONIER, such a detector would 
increase the sensitivity by 1.5 mag and would solve the read- 
out speed and saturation issues. 
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- PIONIER has the status of a visitor-instrument at Paranal 
Observatory. The possibility to offer PIONIER to the com- 
munity is currently investigated. 
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Appendix A: Observations of 8 Sco 

Table A.l. Log of the observations of 6 Sco. The number of observations is reported taking into account the spectral dispersion (3 
spectral channels). Each measurement (pointing) is generally the average of 3 or 5 files of 100 scans each. 
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Fig. A.l. Observed calibrated visibilities (V 2 ) on 6 Sco versus the spatial frequency (m/yum). Each panel is for a baseline. The red 
curves are the result of the global fit by a disk+point model. The (w, v) plane coverage (m/yum) of the observations is displayed in the 
bottom-right panel. North is up and East is to the right. The radial tracks come from the spectral dispersion over 3 spectral channels. 
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Fig. A.2. Observed calibrated closure phase (deg) on S Sco versus the largest spatial frequency of the baseline triple (m/yum). Each 
panel is for a baseline triplet. The red curves are the result of the global fit by a disk+point model. 
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Appendix B: Observations of HIP11231 

Table B.l. Log of the observations of HIP1 1231. The number of observations is reported taking into account the spectral dispersion 
(7 spectral channels). Each measurement (pointing) is generally the average of 3 or 5 files of 100 scans each. 



50 



-50 



_ 1 1 















i 






\ ♦/ 




_ . i 


1 . " 



Date 


Average MJD 


Baseline 


#of V 2 


# closure 


Time spend 


2010-10-29 


55499.2 


D0-H0-G1-I1 


24x7 


16x7 


2.2h 


2010-10-30 


55500.2 


D0-H0-G1-I1 


18x7 


12x7 


1.4h 


2010-10-31 


55501.1 


D0-H0-G1-I1 


18x7 


12x7 


1.2h 


2010-11-30 


55531.1 


E0-G0-H0-I1 


12x7 


8x7 


lh 


2010-12-01 


55532.1 


E0-G0-H0-I1 


12x7 


8x7 


0.8h 


2010-12-02 


55533.0 


E0-G0-H0-I1 


12x7 


8x7 


0.9h 


2010-12-20 


55551.1 


A0-K0-G1-I1 


12x7 


8x7 


0.65h 


-29 


2010-10- 


-30 


2010- 


-10-31 





-50 



50 -50 



50 -50 



2010-12-01 



2010-12-02 



2010-12-20 



50 



-50 



_ 1 1 


i 1 _ 


_ 1 i 




<s> * 






















* 


i 


i r 




i 



1 i 


i 1 _ 

/ 


/ 

/ 

i 


i ,: 



-50 



50 -50 



50 -50 



50 



2010-11-30 



50 -50 



50 



Fig. B.l. (u, v) plane coverage (m/jum) of the observations of HIP1 1231. North is up and East is to the right. The radial tracks come 
from the spectral dispersion over 7 spectral channels. Each panel is for an epoch. 
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Fig.B.2. Observed calibrated visibilities (V 2 ) on HIP11231 versus the spatial frequency (m/yum). The red curves are the result of 
the global fit by a binary model. Each panel is for an epoch. 
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Fig.B.3. Observed calibrated closure phase (deg) on HIP11231 versus the largest spatial frequency of the baseline triple (m/yum). 
The red curves are the result of the global fit by a binary model. Each panel is for an epoch. 
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